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ABSTRACT 

A study was conducted at the NASA Lewis Research Center to obtain a correlating 
parameter which would relate the rate of pressure rise to the volume of spherical liquid- 
hydrogen tankage. A nonventing 22 -inch (56 -cm) diameter spherical tank partially filled 
with liquid hydrogen was subjected to four uniformly heated self -pressurization tests 
with various combinations of percent filling and heat-transfer rate. These data were 
compared with data from a similar study on a 9 -inch (23 -cm) diameter, uniformly heated, 
spherical tank. The two sets of data verified the analysis which predicted the same 
change in pressure in each tank for equal values of heat added per unit volume. 
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EFFECT OF SIZE ON NORMAL-GRAVITY SELF-PRESSURIZATION 
OF SPHERICAL LIQUID HYDROGEN TANKAGE 
by John C. Aydelott and Charles M. Spuckler 
Lewis Research Center 

SUMMARY 

A study was conducted at the NASA Lewis Research Center to obtain a correlating 
parameter which would relate the rate of pressure rise to the volume of spherical liquid- 
hydrogen tankage. A nonventing 22 -inch (56 -cm) diameter spherical tank partially filled 
with liquid hydrogen was subjected to four uniformly heated, self -pressurization tests 
with various combinations of percent filling and heat -transfer rate. These data were 
compared with data from a similar study on a 9-inch (23 -cm) diameter, uniformly heated, 
spherical tank. The two sets of data verified the analysis which predicted the same 
change in pressure in each tank for equal values of heat added per unit volume. 


INTRODUCTION 

Hydrogen is very attractive for use as a rocket fuel because of the high specific im- 
pulse it produces either when reacted with an oxidizer or when used as a heated expel - 
lant as in a nuclear rocket. However, mission planners and tank designers are confronted 
with problems that result because hydrogen has many properties which set it apart from 
common fluids. Gaseous hydrogen has a low density, so that it is practical, for space 
missions, to store hydrogen only as a liquid or solid. The liquid and solid, however, have 
a very low equilibrium temperature, and this makes it very difficult to insulate a storage 
container sufficiently to prevent a net heat gain. As heat is continually added to the tank, 
the tank pressure will rise until venting is required, with an accompanying propellant loss. 

The rate of pressure rise in a hydrogen storage tank is highly dependent on the ther- 
mal transport properties of hydrogen, which make it possible for subcooled liquid and 
highly superheated vapor to coexist in the same container. The temperature of the inter- 
face between the liquid and vapor phases increases because the interface temperature is 



always equal to the saturation temperature corresponding to the increasing total system 
pressure. The average liquid temperature may increase at a slower rate, thus the liquid 
bulk becomes subcooled; the average vapor temperature may increase at a faster rate, 
thus the vapor becomes superheated. The rate of liquid subcooling and vapor superheating 
will vary with the experimental conditions (tank geometry and percent filling, heat- 
transfer rate and distribution, etc. ) and cause the rate of pressure rise to vary. Con- 
sequently, simple thermodynamic analysis often cannot predict the rate of pressure rise 
in a closed system containing liquid hydrogen, and the possibility of developing scaling 
parameters should be considered. 

The majority of the work in this field has been restricted to cylindrical tanks with 
heating only on the side walls. Therefore, natural convection theory for vertical plates 
could be used to predict the heat and mass transfer within the liquid phase. In general, 
direct heating of the vapor was not considered. An excellent review of this work can be 
found in reference 1. Of more direct application to understanding the energy transfer 
processes which take place in any liquid -hydrogen container is the effect of bottom heat- 
ing (refs. 2 to 5). As the proportion of bottom heating increases, turbulent convection 
and boiling cause the temperature of the liquid to become more uniform. 

In contrast to the available information on cylindrical tanks, very little experimental 
or theoretical information is available on the normal -gravity self -pressurization of 
spherical liquid -hydrogen tankage. The earlier work of the first author (ref. 6) deter- 
mined the importance of the heating configuration on the thermodynamic history of spher- 
ical, 9 -inch (23 -cm) diameter, liquid -hydrogen tankage. Heating the tank only from the 
top caused pressure rise rates approximately four times greater than those which resulted 
from heating only the bottom of the tank. Reference 7 presents the results of a self- 
pressurization test on a 50 000-gallon (189 000-liter) spherical liquid-hydrogen storage 

tank subjected to a uniform heat-transfer rate of approximately 0. 6 Btu per hour per 

o 

square foot (1. 9 W/m ). The measured rate of pressure rise was approximately 10 times 
greater than that which would result if the contents of the tank were homogeneous (uniform 
temperature throughout). The lack of experimental data and the difficulty in analytically 
relating the rate of pressure rise to the energy input for full-size hydrogen tanks suggest 
the use of a scaling parameter. 

This study was conducted to obtain a correlating parameter for the effect of size on 
normal -gravity self -pressurization of spherical liquid -hydrogen tankage. This report 
presents the information obtained from four self -pressurization tests of a spherical 
22 -inch (56-cm) diameter, liquid -hydrogen tank. The tests were conducted under normal- 
gravity conditions with approximately uniform heat addition. The effects of the variables, 
liquid filling (approximately 30 to 80 percent by volume) and average heat -transfer rate 
(approximately 17 to 64 Btu/(hr)(ft^) or 53 to 202 W/m^), on the rate of pressure rise 
were determined. The data are compared with theoretical analyses and with the results 
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of reference 6 (9 -in. (23 -cm) diam tank) to show the effect of size on the self- 
pressurization of spherical liquid -hydrogen tankage. The information presented in ref- 
erence 7 (50 000 gal (189 000 liter) tank) is discussed in light of the results of the two 
smaller tank studies. 


ANALYSIS 

Thermodynamic Pressure Rise 

The first law of thermodynamics is 


Q= AU+/p dV 

(Symbols are defined in appendix A. ) For a closed, nonexpanding system, all the heat 
absorbed by the system manifests itself in a change in the total internal energy of the 
system; that is, for dV = 0 


Q = AU 

If the system in question is a tank containing a liquid and its vapor, a knowledge of how 
the added heat affects the internal energy distribution and, thus, the temperature dis- 
tribution within the tank permits prediction of the total system pressure. This is possi- 
ble because, for a liquid-vapor system, temperature and pressure are interrelated at 
the interface (ref. 1). 

Since the temperature distribution in the fluid within a cryogenic storage tank is 
highly complex and is affected by many variables, the most important of which are tank 
geometry, percent filling, heat -flux rate, and heat -flux distribution, any analytical pre- 
diction of pressure rise as a function of heat added must be based on simplifying assump- 
tions. This report makes use of two simple pressure rise models. The mathematical 
development of each of the models is given in appendix B. The first model assumes 
homogeneous conditions (uniform temperature) throughout the hydrogen container and is 
a common calculation that is performed to compare data of this type. The second model 
assumes that all the energy absorbed by the hydrogen container goes into evaporating 
liquid and maintaining the vapor at the saturation temperature corresponding to total - 
system pressure. The liquid -phase temperature remains constant at the saturation tem- 
perature corresponding to the initial pressure. These models are not intended to be 
attempts at describing the process that actually takes place in a nonventing hydrogen con- 
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tainer, but are intended to be a means of comparing one set of experimental data with 

another. The position that experimental data assume in relation to the theoretical models 

on a plot of total pressure against heat added is then a qualitative indication of how the 

energy is being distributed within the hydrogen container. 

Figure 1 is a theoretical plot of tank pressure as a function of heat added for the two 

3 

energy distribution models. The plot is for a 1 -cubic -foot (0. 02 832 -m ) container and 
initial liquid fillings of 25, 50, and 75 percent by volume. 


Similarity Relation 

Starting with the previously presented first law of thermodynamics for a closed non- 
expanding system 


Q = AU = A (mu) = A(pu)V 

— = A(pu) 

V 

For any fluid the product of density and specific internal energy can be expressed as a 
function of the pressure and temperature. However, if the fluid is saturated, the pres- 
sure and temperature are dependent variables, and the vapor mass fraction of the fluid is 
also required in order to specify the desired properties 

pu = f(P, T, X) 

— = — AP + — AT + — AX 
V 9P 9T ax 

If we wish to compare two systems containing only liquid and gaseous hydrogen, 0f/0P, 
0f/0T, and 0f/0X will be the same for each system. Consideration will be limited to 
systems of identical geometry, initial percent filling, and heat -transfer rate in order to 
logically expect similar processes to occur in the two systems. From the last equation, 
it can be seen that similar changes in pressure, temperature, and vapor mass fraction 
may then be expected in each tank for equal values of heat added per unit volume. Bailey 
(ref. 8) reached an identical conclusion. 

This similarity relation holds exactly for both of the models presented in the Thermo- 
dynamic Pressure Rise section. The reader may approximate the energy input, as de- 
termined by these models, that will cause a specified change in pressure for any tank 
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size or filling. Referring to figure 1, this is done by interpolating to determine the ef- 
fect of percent filling and by multiplying the heat added by the volume of the tank in cubic 
feet. 


APPARATUS 

Liquid Hydrogen Container 

Figure 2 is a cross-sectional drawing of the liquid -hydrogen apparatus that consisted 
of three concentric spheres; the inner sphere contained the liquid hydrogen, the interme- 
diate sphere had electric heating coils mounted on its exterior surface, and the outer 
sphere served as a vacuum jacket to reduce the gaseous conduction of heat. The outer 
surface of the inner sphere and the inner surface of the heaters were painted black in 
order to increase their emissivity. The stainless -steel fill and vent tubes supported the 
inner sphere. 

A heater controller, which basically consisted of a bridge circuit which balanced the 
resistance of a temperature sensor on each heater with a corresponding rheostat on the 
control panel, was used to maintain heater temperatures of 360° or 500° R (200 or 
278 K). 


Instrumentation 

Temperature and pressure transducers measured the total system pressure, vacuum - 
space pressure, surface temperature of the inner sphere, heater, and vacuum jacket, and 
temperature at 16 locations inside the inner sphere. Figure 3 shows the location of the 
platinum temperature transducers on the inner sphere and the four carbon resistor 
temperature rakes that were located within the inner sphere to measure the temperature 
of the hydrogen liquid and vapor. Figure 4 shows two of the carbon resistor rakes. 

At any time prior to a test run, the resistance of any temperature transducer could 
be determined by the use of a digital ohmeter mounted in the control panel. A very small 
electric current was passed through a temperature transducer to determine the resistance 
and thus the temperature of the transducer. When a relatively high current was applied 
to a carbon resistor, its temperature was quite different depending on whether the tem- 
perature probe was in the liquid or vapor phase. This was due to self -heating. Exploita- 
tion of this fact, together with careful arrangement of the carbon resistors, made it pos- 
sible to use the carbon resistors to determine the liquid level in the sphere prior to the 
start of a test. 
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A detailed description of the temperature and pressure measuring and recording sys- 
tems and an analysis of the error associated with them will be found in reference 6. 


PROCEDURE 

Prior to the assembly of the experiment, thermocouples were attached to the inner 
sphere, heaters, and vacuum jacket. All temperature transducers were calibrated by 
submerging the three spheres in a liquid nitrogen bath at 139. 5° R (78 K) and by placing 
the spheres in a carefully controlled temperature chamber at 347°, 486°, 537°, and 601° R 
(193, 270, 298, and 334 K). After the experiment had been assembled, the inner-sphere 
temperature transducers were calibrated at 36. 4° R (19. 1 K) by filling the sphere with 
liquid hydrogen. Each bridge was calibrated by using a decade box to obtain a voltage 
versus resistance plot. Prior to each test, the pressure transducers were calibrated 
with standard pressure gages. 

Each of the calibration curves for the temperature transducers and for the bridges 
was curve fitted using a digital computer. The magnetic data tapes from each test were 
fed into the digital computer along with the calibration curve fits, and an automatic data 
reduction program returned printed temperature data at half -second intervals for each 
transducer. A continuous plot of the output of the pressure transducers was obtained by 
feeding the output of the magnetic data tapes into a line recorder. 

For each of the tests, the experiment was prepared in an identical manner: only the 
actual test conditions were varied. The space between the inner and outer spheres was 
evacuated first with a mechanical pump and then with a diffusion pump. Then the inner 
sphere was filled with liquid hydrogen. The addition of the liquid hydrogen reduced the 
pressure in the space between the inner and outer spheres due to cryogenic pumping. The 
capillary tubes (fig. 2) were used to remove liquid from the inner sphere until the approx- 
imate desired liquid level was obtained. 

For the low heater temperature (low heat-transfer rate) tests, a gas meter installed in 
the vent line together with the carbon resistors in the inner sphere (described in the In- 
strumentation section) made it possible to determine the liquid level at the beginning of the 
test. The resistors accurately determined the liquid position at some time prior to the 
beginning of the test. At this time the gas meter began to record the volume of hydrogen 
vapor which then leaves the inner sphere. Measurements of the vapor temperature at the 
gas meter and the atmospheric pressure determine the density of the vapor. From a know- 
ledge of the volume and density of the vapor, the mass of hydrogen that leaves the inner 
sphere before the test begins is calculated and the initial percent filling can be determined. 

For the high heater temperature (high heat -transfer rate) test, the boiloff rate of 
the hydrogen was greater than the capacity of the gas meter. Consequently, the eight 
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carbon resistor temperature transducers that were located near the 50-percent fi llin g; 
level were used to determine the initial filling. By recording the time when the liquid 
level went below each resistor and extrapolating to the starting time of the test, the ini- 
tial percent filling could be determined. 

The heater controller was set to maintain the desired heater temperature prior to 
determining the liquid level to eliminate transients. At minus 1 minute, the system be- 
gan recording data on magnetic tape. At zero time the vent valve was closed and the ex- 
periment was allowed to self -pressurize. The tests were terminated either at a nomi- 

o 

nal pressure of 100 psia (68. 95 N/cm ) or after 45 minutes if the pressure did not reach 
100 psia (68. 95 N/cm\ After the experiment reached 100 psia (68.95 N/cm^) or ran 
for 45 minutes, the vent valve was opened and the pressure was allowed to decay slowly. 
If a sufficient amount of liquid hydrogen still remained in the experiment, a similar test 
at a lower filling was run as soon as the new liquid level had been determined. 


EXPERIMENTAL RESULTS 

Four quiescent tests with uniform heating were performed. Nominal initial fillings 
of 30, 50, and 80 percent for the low heat -transfer -rate tests and 50 percent for the high 
heat -transfer -rate test were desired. Heater temperatures of approximately 360° and 
500° R (200 and 278 K) were used for the tests. 

Figures 5 to 8, showing pressure and temperature as a function of time, present the 
data obtained from the test runs. Only enough data symbols are included to identify the 
curves. Each of the figures is made up of four plots: (a) total pressure as a function of 
time, (b) outer -sphere and average heater temperature as a function of time, (c) upper 
inner -sphere internal and surface temperature as a function of time, and (d) lower inner - 
sphere internal temperature as a function of time. Figures 5 to 7 show the effects of three 
different percent fillings (31. 6, 48. 9, and 79. 8 percent) for a nearly constant heat -flux, 
approximately 20 Btu per hour per square foot (63 W/m ). Figures 6 and 8 show the ef- 
fects of two different average heat fluxes, 18. 9 and 64. 4 Btu per hour per square foot (60 
2 

and 203 W/m ), for an approximately half -full sphere. 


DISCUSSION OF EXPERIMENTAL RESULTS 
External Heat Transfer 

A first look at the data might suggest that the rate of pressure rise is the most im- 
portant parameter to consider when examining a group of tests which had identical heater 
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temperatures. However, because the hydrogen vapor becomes superheated and because 
the top of the inner sphere increases in temperature as a test proceeds, this approach can 
be quite misleading. The heating of the inner sphere increases the temperature of the 
container wall causing a reduction in the rate of energy input to the hydrogen. This is due 
to both the reduction in radiant exchange and the energy that is absorbed by the container 
wall. Higher heater temperatures and lower percent liquid fillings cause increasing 
inner -sphere temperatures so that the average heat flux to the hydrogen is not a function 
only of heater temperatures. As a result, it is better to compare the amount of heat that 
must be added to the hydrogen container to cause a given pressure rise for a particular 
set of conditions. The method of calculating the net heat addition to the hydrogen container 
used in these experiments is given in appendix C. The results of this analysis indicate 
that the main source of energy input to the experiment was radiant exchange from the 
heaters. Heat transfer due to conduction along the fill and vent tubes and instrumentation 
wires, as well as gaseous conduction through the vacuum space, was small. 

Table I is a summary of the experimental results and the heat-transfer analysis for 
both the present study and the uniform heating tests of the 9 -inch (23 -cm) diameter 
spherical tank (ref. 6). The initial percent filling was determined as explained in the 
PROCEDURE section. The pressure -rise rate is an average value obtained by dividing 
the change in the pressure during the test by the total test time. The bulk temperature 
was assumed to be the lowest recorded temperature from the lower carbon resistor 
temperature transducers. Dividing the change in the bulk temperature by the change in 
the saturation temperature during the test gives an indication of how much energy went 
into heating the liquid and, thus, how nearly homogeneous the liquid is at the end of 
the test. The maximum change in the vapor temperature is an indication of how much 
energy went into superheating the vapor. However, the vapor temperature data from 
the first two tests on the larger spherical tank should not be compared with the other 
tests. If the first two tests on the larger tank had not been terminated before reaching 
a pressure of 100 psia (68. 95 N/cm ), the maximum change in the vapor temperature 
would have been larger. The average heat flux is determined by dividing the total energy 
input to the hydrogen as calculated in appendix C by the test time and the surface area 
of the inner sphere. This average is then broken down into the heat flux through the 
liquid wetted walls and through the walls exposed to vapor by assuming that none of the 
reradiated heat flux (eq. (C8)) comes from the liquid wetted walls. Breaking the heat 
flux into parts in this manner clearly shows the effect of the increasing upper inner-sphere 
temperature which reduces the net radiant heat exchange. 
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Pressure-Rise Characteristics 


Effect of per cent filling . - Figure 9 shows the effect of percent filling on the sphere 
pressure as a function of heat added for six uniformly heated tests. Figure 9(a), taken 
from reference 6 (9 -in. (23 -cm) diam spherical tank), and figure 9(b) for the present 
study (22 -in. (56 -cm) diam spherical tank), show similar results. The pressurization 
rate was only slightly affected by varying the percent filling with a trend toward higher 
pressure -rise rates at higher fillings. In order to understand why the pressure -rise 
rate is only slightly affected by the percent filling, reference is made to table I. The 
data for the tests indicate that, as the percent filling is increased, the liquid becomes 
more subcooled (less homogeneous) and the vapor becomes less superheated (more 
homogeneous). These two effects tend to counterbalance each other. 

Effect of heat -tran sfer rate . - Figure 10 shows the effect of heat -transfer rate on 
sphere pressure as a function of total heat added for the approximately 50 -percent filled, 
uniformly heated, quiescent tests. Figure 10(a), taken from reference 6, and figure 
10(b), for the present study, show similar results. Since the two coordinates, pressure 


TABLE I. - SUMMARY OF EXPERIMENTAL RESULTS AND HEAT -TRANSFER ANALYSIS 3, 


Test 

Initial 

Average 

Liquid - 

Maximum vapor 

Unit area heat -transfer rate q/A- 


filling, 

pressure -rise rate 

temperature 

temperature change 








percent 

— ■ 

/ 2 

ratio, 

o„ 


Average 

Wetted area 

Dry area 



psi 

N/ cin 

A V AT sat 

R 

K 

Btu 

W 

Btu 

W 

Btu 

W 



min 

min 
















(hr) (ft 2 ) 

2 

m 

(hr)(ft 2 ) 

m 2 

(hr)(ft 2 ) 

m 2 




22 -in. 

(56 -cm) diameter spherical tank -uniform heating 




1 

31.6 

1.23 

0. 85 

0. 59 

b (130) 

b (72) 

16. 9 

53 

22. 3 

70 

13.6 

43 

2 

48.9 

1.41 

.97 

.48 

b (155) 

b (86) 

18.9 

60 

22. 0 

69 

15. 8 

50 

3 

79. 8 

2.04 

1.4 

. 34 

164 

91 

21.9 

69 

22.6 

71 

20.2 

64 

4 

54. 2 j 

4. 78 

3.3 

. 51 

222 

123 

64.4 

203 

83. 7 

264 

41.6 

131 




9 -in. (23 -cm) diameter spherical tank -uniform heating (ref. 6) 




1 

34.6 

3. 5 

2.4 

0. 55 

110 

61 

17. 5 

55 

20.0 

63 

15.7 

49 

2 

51.4 

3.6 

2. 5 

.46 

115 

64 

18.2 

57 

20.0 

63 

16.4 

52 

3 

34.9 

11. 3 

7. 8 

.62 

212 

118 

59.9 

189 

80.3 

253 

45. 8 

144 

4 

48. 9 

13. 5 

9.3 

.47 

200 

111 

65.0 

205 

80. 7 

254 

48.9 

154 

5 

76. 5 

17. 0 

11.7 

.32 

159 

88 

72.6 

229 

80.5 

254 

54. 1 

171 

6 

36. 8 

19. 3 

13.3 

.61 

258 

143 

105.4 

332 

148.0 

466 

74. 5 

235 

7 

50. 7 

23.6 

16.3 

.47 

245 

136 

111.6 

352 

142.7 

450 

78.4 

247 

8 

77.2 

30.3 

20.9 

.31 

182 

101 

128.6 

405 

141.9 

447 

97. 7 

308 


a Quantities defined in heat transfer analysis and test results section. 

b 2 

Test terminated before reaching 100 psia (68. 95 N/cm ). 
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and heat added, are the integrals over time of pressure-rise rate and heat -transfer rate, 
coincident test data indicate a linear relation between pressure -rise rate and heat- 
transfer rate; that is doubling the heat-transfer rate will double the pressure-rise rate. 
This linear relation was followed almost exactly for all the approximately 50-percent 
filled, uniformly heated, quiescent tests. 


Internal Heat and Mass Transfer 

Liquid bulk. - Analysis based on the modified Rayleigh number (Ra = C aSp^qD^/ 

2 IQ IQ P 

k fi A) for the present study and the work of reference 6 (Ra » 10 to 10 ), indicates 
that the mode of heat transfer in the liquid bulk would be turbulent convection (refs. 9 
to 11). However, a summary of liquid -hydrogen boiling studies presented in references 
12 and 13 indicate that, at the heat fluxes employed for these tests, nucleate boiling is 
quite likely to occur. An essentially uniform -temperature liquid bulk would be antici- 
pated for either turbulent convection or boiling heat transfer. Because a uniform - 
temperature liquid bulk was experimentally observed, it is probable that, at the lower 
heat fluxes, the heat and mass transfer in the liquid was dominated by turbulent convec- 
tion with some boiling entering in at the higher heat fluxes. 

In contrast to the results obtained with the smaller tanks, is the data (ref. 7) obtained 
from a much larger spherical tank (50 000 gal (189 000 liter)). For this test the heat- 
transfer rate was approximately 0. 6 Btu per hour per square foot (1. 9 W/m ). The 
pressure -rise rate was approximately 10 times greater than the calculation based on 
homogeneous conditions, which is much greater than for any of the small tank tests. It 
is believed that, due to the low heat-transfer rate, a laminar natural convection boundary 
layer along the tank wall carried the incoming energy directly to the liquid surface. This 
conclusion is based on the fact that very little heating of the liquid bulk was observed dur- 
ing the experiment. Natural convection boundary-layer flow would be expected if the Ray- 
leigh number for this test was less than the Rayleigh number for the small tank tests. 
However, the Rayleigh number, based on tank diameter, for the reference 7 test is much 
larger due to the fourth power dependence of the Rayleigh number on the significant di- 
mension. Consequently, either this type of Rayleigh number analysis is not valid or the 
characteristic dimension in the Rayleigh number should be some other quantity such as 
the boundary -layer thickness. 

Vapor. - The heat -transfer processes which take place in the vapor are not clearly 
understood. The fact that the lines of constant temperature in the vapor were horizontal, 
with increasing temperature at higher vertical positions, may rule out the possibility of 
any convective flow. The validity of this point is shown by figure 11, which is a typical 
plot of inner -sphere temperature as a function of vertical position. The liquid -vapor in- 


10 


terface is located at a height to radius ratio of approximately 1. This figure shows that 
for any given time, the temperature of the vapor space was only a function of the vertical 
coordinate, since the data for all the instrumentation, both centrally located and near the 
container wall, had the same temperature profile. It is possible that an involved conduc- 
tion analysis could predict the vapor temperature gradients; however, the results of the 
study presented in reference 14 indicate that the increasing system pressure could cause 
fluid motion which would dominate any conduction effects. 

L iquid-vapor interface . - The transfer of heat and mass across the liquid -vapor 
interface, either by evaporation or condensation, is also a process that is not well under- 
stood. Reference 15 presents an analysis based on irreversible thermodynamics that 
will predict the behavior of a fluid as it changes phase. However, the three transport 
coefficients that are necessary for solution of the equations must be experimentally de- 
termined. To the authors knowledge, the values of these coefficients for hydrogen are 
not available. 

Liquid -thermal layer . - The liquid -thermal layer is the variable temperature 
region in the liquid between the liquid -vapor interface and the uniform temperature liquid 
bulk. Since the liquid bulk was always experimentally observed to be subcooled and the 
liquid-vapor interface is always saturated, it is reasonable to assume that conduction 
heat -transfer could account for the temperature gradients that exist in the liquid thermal 
layer. Here analysis is difficult because of the wide diversity of available data for the 
thermal conductivity of liquid hydrogen (ref. 16). 


Effect of Tank Size 

As can be seen, a great deal more basic research is necessary before the internal 
heat and mass -transfer processes which take place in a closed cryogenic container can 
be predicted analytically. Until this information is available, the rate of pressure rise 
in a closed cryogenic container must be determined experimentally. However, the use 
of models is a common engineering tool that can be used once scaling laws have been 
derived. The ANALYSIS section shows that the same pressure rise should result in two 
tanks of similar geometry, percent filling, and heating configuration if the ratio of the 
total heat added to tank volume is the same. The effect of large differences in heat- 
transfer rate was shown by the comparison of the results of reference 7 (50 000-gal 
(189 000-liter) spherical tank) with the results obtained in the smaller tanks. The large 
tank had a rate of pressure rise approximately 10 times greater than would result from 
the homogeneous analysis while the rate of pressure rise in the smaller tanks was three 
to four times greater than homogeneous. It is believed that this difference in compara- 
tive pressure -rise rates is due to the large difference in heat -transfer rates between the 
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larger and small tanks. Consequently, for scaling purposes, the heat -transfer rates 
must be such that similar heat- and mass-transfer processes take place in the liquid 
bulk. Figure 12 is a plot of pressure as a function of total heat added divided by tank 
volume for two pairs of similar tests performed using the two smaller spherical liquid- 
hydrogen tanks. These two pairs of data were chosen because their initial fillings were 
nearly identical, thus fulfilling one of the restrictions postulated in the ANALYSIS sec- 
tion. These data fell within a narrow band. However, as shown in figure 9, the effect of 
percent filling was small, so that all the data would fall within a narrow band, thus sup- 
porting the analysis which predicted that the effect of size on the self -pressurization of 
hydrogen tankage can be expressed as a simple geometric relation. 


Energy Distribution 


In order to better understand the experimental results, an analysis was performed 
to determine what percentage of the incoming energy resulted in heating of the liquid 
bulk, the liquid-thermal layer, evaporation of liquid, and superheating of the vapor. The 
details of these energy -distribution calculations will be found in appendix D. Table II is 
a summary of the results of the energy distribution analysis for both the present study 


TABLE E. - RESULTS OF ENERGY -DISTRIBUTION ANALYSIS 


Test 

Average 

Energy input to 

Liquid 

Vapor 

Evapo- 


filling, 

percent 

liquid wetted walls, 
percent of total 

Total 

Bulk 

Thermal 

layer 


ration 




Heat added, percent of total energy input 


22 -in. (56 -cm) diameter spherical tank 

-uniform heating 


1 

32.6 

50.0 

75. 8 

45.4 

30.4 

15.0 

9.2 

2 

50. 1 

58.3 

81.9 

64.7 

17.2 

11.4 

6.7 

3 

81. 5 

74.2 

93.9 

82.7 

11.2 

4.8 

1.3 

4 

56.0 

70. 5 

82.9 

67.6 

15.3 

14.0 

3. 1 

9 -in. (23 -cm) diameter spherical tank -uniform heating (ref. 6) 

1 

35. 8 

46.4 

77.4 

32. 5 

44.9 

8. 8 

13.8 

2 

52.9 

57.2 

83.0 

51.7 

31.3 

6.7 

10.3 

3 

36.7 

54.9 

76. 1 

39. 1 

37.0 

18. 1 

5. 8 

4 

50.4 

62.7 

79.4 

58.4 

21.0 

15.4 

5.2 

5 

78.2 

77. 5 

91. 1 

64.6 

26. 5 

7.7 

1.2 

6 

38.4 

59.2 

74.4 

47. 8 

26.6 

22.0 

3.6 

7 

52.2 

66.0 

77.4 

65.5 

11.9 

18.7 

3.9 

8 

78.7 

77. 1 

90.7 

63.8 

26.9 

8.7 

.6 
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and the uniform -heating tests of reference 6. The analysis shows that the energy input 
to the dry walls heated the vapor, supplied the necessary energy for evaporation, and 
heated some of the liquid since the energy input to the liquid wetted walls was always less 
than the input to the liquid. The energy input to the liquid wetted walls heated the liquid 
bulk, and in some cases, a portion of the liquid -thermal layer. Essentially, the same 
conclusion was reached by the authors of reference 17 as a result of their tests on a 
625-gallon (2360-liter) cylindrical hydrogen Dewar. 

The pressure -rise rate in the hydrogen container for the uniformly heated tests is 
primarily a function of the energy input to the dry walls since the energy input to the 
liquid wetted walls only heats the liquid. The only contribution that the heating of the 
liquid wetted walls makes to the container pressure is due to the thermal expansion of 
the liquid. Liquid hydrogen does have a relatively high coefficient of thermal expansion, 
but this effect is secondary to the energy input to the vapor and for evaporation in deter- 
mining the rate of pressure rise. 

The minor differences in the energy input distribution explain the slight deviations 
in the results that were obtained from the two small hydrogen tanks (fig. 12). As the 
heat-transfer rate was increased, the contents of the 9 -inch (23 -cm) tank became slightly 
less homogeneous (fig. 10(a)); while the contents of the 22 -inch (56 -cm) tank remained 
essentially the same (fig. 10(b)). This difference is attributed to the thicker wall of the 
22 -inch (56 -cm) tank, which absorbed a larger proportion of the incoming energy at the 
highest heat -transfer rates and caused less heating of the vapor (see tables I and II). 


SUMMARY OF RESULTS 

A nonventing 22 -inch (56 -cm) diameter spherical tank partially filled with liquid 

hydrogen was subjected to four uniformly heated self -pressurization tests. These tests 

were conducted with various combinations of liquid filling (approximately 30 to 80 percent 

o 

by volume) and heat -transfer rate (approximately 17 to 64 Btu/(hr)(ft ) or 53 to 202 
W/m ). These data were compared with data from a similar study performed on a 
9 -inch (23 -cm) diameter, uniformly heated, spherical tank (ref. 6) and with the self- 
pressurization of a 50 000-gallon (189 000-liter) spherical liquid-hydrogen storage tank 
(ref. 7). For the range of variables investigated, the following results were obtained: 

1. The data from the 9 -inch (23 -cm) and 22 -inch (56 -cm) diameter spherical tank 
tests verified the analysis that predicted that the effect of size on self -pressurization of 
hydrogen tankage can be expressed as a simple geometric relation involving the heat 
added per unit volume. This conclusion is supported by the fact that the pressure data, 
plotted as a function of heat added per unit volume, fell within a narrow band. 
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2. The pressure -rise rate increased almost linearly with increasing heat -transfer 
rate; however, as the heat-transfer rate was increased, the contents of the 9 -inch (23-cm) 
tank became slightly less homogeneous while the contents of the 22 -inch (56 -cm) tank 
remained essentially the same. This difference is attributed to the thicker wall of the 

22 -inch (56 -cm) tank which absorbed a larger proportion of the incoming energy at the 
higher heat -transfer rates and caused a change to more bottom heating. 

3. For both the 9 -inch (23-cm) and 22 -inch (56 -cm) tanks the pressure rise rate was 
only slightly affected by varying the percent filling with a trend toward increasing 
pressure -rise rates at higher fillings. 

4. The 50 000-gallon (189 000-liter) spherical liquid -hydrogen storage tank (ref. 7) 
had a much greater rate of pressure rise, when compared with the homogeneous analysis, 
than either of the smaller tanks. This difference was attributed to a different mode of 
heat transfer in the liquid. At the low heat fluxes encountered in the reference 7 experi- 
ment (0. 6 Btu/(hr)(ft ) or 1. 9 W/m ), heat is transported in the liquid by laminar nat- 
ural convection rather than turbulent convection as is the case for the tests on the smaller 
tanks so that larger temperature gradients and nonhomogeneous conditions result. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio December 12, 1968, 

124-09-17-01-22. 
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APPENDIX A 


SYMBOLS 


2 2 

A surface area, ft ; m 

a accommodation coefficient 

absorption factor 

C specific heat, Btu/(lb)(°R); 

J/(kg)(K) 

Cy specific heat, Btu/(lb)(°R); 

J/(kg)(K) 

D diameter, ft; m 

F angle factor 

g acceleration due to gravity, 

ft/sec 2 ; m/sec 2 

k thermal conductivity, Btu/ 

(hr)(ft)(°R); W/(m)(K) 

L length, ft; m 

M molecular weight 

m mass, lb; kg 

2 

P pressure, psia; N/cm 

pf filling by volume, percent 

Q heat added, Btu; J 

q heat -transfer rate, Btu/hr; W 

R universal gas constant, 

1.986 Btu/(lb)(mole)(°R); 

8. 3143 J/(mole)(K) 

r reflectivity 

T absolute temperature, °R; K 

t time, hr 

U total internal energy, Btu; J 


u specific internal energy, Btu/Ib; 

J/kg 

V volume, ft**; m 8 

X vapor mass fraction 

x linear distance, ft; m 

/3 coefficient of thermal expansion, 

l/°R; 1/K 

y ratio of specific heats, C /C y 

P 

e emissivity 

o 

p dynamic viscosity, (lb)(sec)/ft ; 

(N)(sec)/m 2 

p density, lb/ft 8 ; kg/m 8 

a Stefan -Boltzmann constant, 

0. 1713X10" 8 Btu/ (hr)(ft 2 )(°R 4 ); 
5. 6697X10' 8 W/(m 2 )(K 4 ) 

Subscripts: 

a absorbed 

av average 

b liquid bulk 

ev evaporation 

f final or any intermediate state 

gc gaseous conduction 

i initial state 

j summation variable 

L at x = L 

L liquid 

l liquid thermal layer 
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m 

mean or average value 

V 

vapor 

max 

maximum 

V g 

at vacuum gage 

n 

summation variable 

w 

wall 

r 

radiant 

0 

at x = 0 

rr 

reradiated 

1 

inner sphere 

s 

system 

2 

upper heater 

sat 

saturation 

3 

lower heater 

sc 

solid conduction 

4 

opening in upper heater 

st 

stored 

5 

combined upper and lower heaters 
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APPENDIX B 


THERMODYNAMIC ANALYSIS 


One method of analyzing a system thermodynamically is to define the conditions at 
the beginning and the end of a process; then the necessary input to the system can be de- 
termined. For the problem of a liquid -hydrogen storage tank, which is sealed at the 
beginning of the test, the initial condition is that the tank contains a homogeneous mixture 
at atmospheric pressure with a known percent filling by volume. The final condition is 
defined by the model being considered. For a nonexpanding closed system, the input is 
heat, and, as stated by the first law of thermodynamics 

Q = AU (Bl) 


Equation (Bl) may be written in the form 


Q = U f - Uj = <m M u M + m y> f u y> f ) - (m^ . + m y> ^ t ) 


(B2) 


The density and specific internal energy of each phase at state i can be found if the sys- 
tem is known to be homogeneous and at the saturation temperature corresponding to at- 
mospheric pressure (ref. 18). The total internal energy at state i can then be deter- 
mined because 


m 




(pf)i 

i V 

100 


(B3) 


v 


m 


v. 1 


= P 





100 


(B4) 


For a closed nonexpanding system, the mass of the liquid plus the mass of the vapor is a 
constant; consequently, the system density is a constant 


(pf)i 

P s = P 

s 100 




l - 


(pOj 


100 


v,i 


(B 5) 


Homogeneous model . - State f (and thus the density and internal energy of each 
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phase) is defined by the fact that the system is homogeneous and is at the saturation tem- 
perature corresponding to the system pressure. Equation (B5), written for state f, can 
be solved for the percent filling at state f 


(Pi) f = (~ ~ Pv ’ f V 00 

W,f p v,f / 


(B6) 


The total internal energy at state f can then be determined using equation (B6) 


and 


m 


l,i~ p £,f 


(pf) f 

100 


(B7) 


m 


V,f 




(B8) 


The amount of heat required to reach state f for the homogeneous model can now be cal- 
culated by using equations (B2) to (B4) and (B7) and (B8). 

Surface -evaporation model. - State i is the same as that for the previous model so 
that the total internal energy at state i is found by the identical procedure. The surface- 
evaporation model is based on the concept that all the energy goes into evaporating the 
liquid and maintaining the vapor at the saturation temperature. If it is assumed that 
the liquid is incompressible, then the density and internal energy of the remaining liquid 
will be unaltered by the process; that is, p £ i = p £ f and u^ . = u £ f . The density and 
internal energy of the vapor are defined by the fact that the vapor is homogeneous and 
at the saturation temperature corresponding to the final system pressure. Equation (B5), 
written for state f, can be solved for the percent filling at state f: 


(pf) f = ff * ~ P v,A 100 


(B9) 


The mass of liquid and vapor at state f can then be determined using equation (B9) 


m 


t,i~ p SL,i 


(pf) f 

Too 


V 


(BIO) 
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and 


m 


r,f- p v,f [* ' 


(p£) f 
100 J 


(Bll) 


The total internal energy at state f can now be determined, and from equations (B 2 ) to 
(B 4 ) and (BIO) and (Bll) the heat required to reach state f for the surface -evaporation 
model can be calculated. 
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APPENDIX C 


HEAT-TRANSFER ANALYSIS 

The amount of energy absorbed by the contained hydrogen is equal to the heat trans - 
ferred to the sphere by radiation, solid conduction, and gaseous conduction minus the 
amount of energy stored in the container itself; that is, 

Q a = <«r + «sc + %c )4t - < C1 > 

The amount of heat transferred by thermal radiation from the heated intermediate 
sphere to the inner sphere is determined by the method presented in reference 19. For 
the radiant exchange calculations the outside of the inner sphere is assigned the number 1, 
the inside of the upper heater, number 2, the inside of the lower heater, number 3, and 
the opening in the upper heater, number 4. The net rate of radiant heat absorbed by the 
inner sphere when both heaters are installed is 

4 

q i = E “ViV? - (c2) 

3=1 

where B^, the absorption factor, is defined as the fraction of the total radiant -energy 
emission of surface j which is absorbed by surface 1. The absorption factors are de- 


termined by solution of the following simultaneous equations: 

( F H r l - l)B n + F i2 r 2 B 21 + F 13 r 3 B 31 + F 14 r 4 B 41 + F ll e l = 0 ( C3 ) 

F 21 r l B ll + ( F 22 r 2 “ ^ B 21 + F 23 r 3 B 31 + F 24 r 4 B 41 + F 21 e l = 0 ( C4 ) 

F 31 r l B ll + F 32 r 2 B 21 + (F 33 r 3 “ 1)B 31 + F 34 r 4 B 41 + F 31 e l = 0 (C5) 

F 41 r l B ll + F 42 r 2 B 21 + F 43 4 3 B 31 + (F 44 r 4 " 1 ) B 41 + F 41 e l = 0 ( C6 ) 


This technique treats all diffuse -radiation circumstances and requires only a know- 
ledge of the geometry of the four surfaces, the average temperature of the surfaces, and 
the emissivity of the surfaces. Lewis Research Center personnel experimentally obtained 
the emissivity of the inner surface of the heaters and the outer surface of the hydrogen 
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container (fig. 13) by using a sample identical to both surfaces. Because of the wires and 
tubes coming through the hole in the upper heater, the heat radiated through the hole from 
the outer sphere was assumed to be emitted from a black body. 

Because of the variation in temperature of the inner sphere, the last term in equa- 
tion (C2) was expressed as an integral and called the inner -sphere reradiated heat flux 

’l, rr = / « T4dA < C7 > 

where is the surface area of the inner sphere. This integral was approximated by 
the summation 


n 

1l,rr“ E «i T J AA i <C8) 

j=l 

It was assumed that the hydrogen temperature profiles were symmetric with respect to 
the vertical axis. In other words, at any time during a test all vertical planes passing 
through the center of the inner sphere would exhibit identical temperature patterns, and 
the left side of such a plane would be the mirror image of the right side. This assumption 
is based on the fact that the sphere, heaters, fill and vent tubes, instrumentation wires, 
and liquid-vapor interface all have symmetry with respect to the vertical axis; conse- 
quently, there is no reason to anticipate that the hydrogen temperature profiles would be 
different on opposite sides of the container. Based on this assumption, horizontal uniform 
temperature sections could be used to divide the sphere into elemental surface areas. A 
digital computer was used to curve fit the emissivity as a function of temperature curve 
and the inner -sphere temperature as a function of position for each time interval. An 
average temperature for each elemental area was used to determine the local emissivity, 
and the summation was performed every 30 seconds by using a digital computer. 

The energy stored in the container at any time interval was also expressed as an 
integral 


Q st = 



pCTdV 


(C9) 


where V. is the volume of the container wall. This integral was approximated by the 
X, w 

summation 
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(CIO) 


n 

Q x « pc -T-AV. 

^st ^ p p,J ] ] 

j=l 

Assuming the same temperature distribution, a digital computer was used to curve fit the 
stainless -steel specific heat as a function of temperature curve found in figure 14 (ref. 20) 
and the inner -sphere temperature as a function of position curves for each time interval. 
An average temperature for each elemental volume was used to determine the specific 
heat and the summation was performed every 30 seconds by using a digital computer. 

The differential equation and boundary conditions for one -dimensional heat transfer 
by solid conduction are the following: 

dx \ dxy 


{ at x = 0, T = Tq 

(Cll) 

at x = L, T = T l 


At the very low temperatures encountered with the use of liquid hydrogen, the thermal 
conductivity of most materials is highly temperature dependent and can be expressed as 
some function of the absolute temperature. Substituting the boundary conditions in equa- 
tion (Cll) and integrating (ref. 19) result in 


^sc 

A 


= k 


m 



(C12) 


where 


k 


m 




k(T)dT 


(C13) 

> 


Figure 15 shows the stainless -steel thermal conductivity as a function of temperature 
(ref. 20) and the curve fit which was used to perform the necessary integration in equa- 
tion (C13). The amount of heat transferred through the vacuum space by gaseous con- 
duction was determined by the following equation (ref. 20). 


22 


i 


■ ■■ mm n ai ii i muni 


I I III III 


III 


III 


J 



q gc _ 1 a l a 5 ? (Tg - T i) (C14) 

A 1 2 a 5 + a l^ " a 5^ A l //A 5 W ” V » 2f ^TygM 
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APPENDIX D 


ENERGY DISTRIBUTION ANALYSIS 

The total volume of the hydrogen container as it appears in this analysis included the 
volume of the plumbing up to the seal -off point of the experiment. Expansion of the con- 
tainer due to the increasing pressure was neglected. 

The total test time was divided into 30-second intervals, and by using an iterative 
method the following calculations were performed for each time interval. As a first ap- 
proximation, it was assumed that the liquid level did not change between successive time 
intervals so that the volume of vapor could be determined. 

Referring to figure 11, the energy distribution analysis was based on the assumption 
that the vapor space could be divided into horizontal uniform temperature disks. Each 
disk or elemental volume approximated a region of constant temperature and pressure. 
Since two thermodynamic properties are known, for each elemental volume, any other 
property can be determined. The properties of particular interest are the density and 
specific internal energy. The equations that were used to generate the hydrogen tables 
(ref. 18) were curved fitted using a digital computer so that the desired properties were 
readily available once the temperature and pressure were obtained from the experimental 
data. By using the elemental volumes to perform a summation (approximating an integra- 
tion) it is possible to determine the total mass and internal energy of the vapor space at 
any time. 


m 


v 



dm 


v 


U = f u dm 
V J V V 


These two integrals were approximated by the summations 


n 


m » V' p . AV ■ 
v H v, j V v, j 

j=l 


(Dl) 


(D2) 


(D3) 
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I 


u * 

v 


n 

E 

1=1 


U 


j p v, 


• AV . 

3 V >J 


(D4) 


At any time these summations can be evaluated by using a digital computer. It is first 
necessary to curve fit the vapor temperature against position data for the test times of 
interest. The computer can then determine the necessary properties, from the hydrogen 
property curve fits, for each elemental volume and perform the necessary mathematical 
operations. The energy input to the vapor during any time interval is 

^v v, f v, 1 


The change in the mass of the vapor is 


Am v =m v ,f -m v ,i 


The energy input that results in evaporation is 


Q gv = (Am v )(Heat of vaporization) 


(D5) 


(D6) 


The energy input to the liquid is determined by subtracting the energy input to the vapor 
and the energy input that results in evaporation from the total energy input to the contained 
hydrogen. 


Q t = Q a - Q v - Q ev <D7) 

To explore further the energy distribution within the liquid phase, the energy input 
to the liquid was broken down into two parts: the energy that went into the bulk of the 
liquid and the energy that went into heating the thermal layer between the saturated 
liquid-vapor interface and the bulk of the liquid. Figure 11 shows that the bulk tempera- 
ture, or lowest recorded temperature, is representative of a large portion of the liquid 
mass. For the purpose of mathematical computation, a linear temperature gradient from 
the bulk temperature to the saturation temperature at the interface was assumed. It is 
realized that for some of the tests this is a poor approximation to the actual temperature 
gradient, but the analysis based on this assumption helps to explain further how energy is 
transported and distributed within the liquid hydrogen. The bulk temperature, at any time 
is determined from the instrumentation and the thermal layer average temperature is 


25 



T = T b + T sat (D8) 

1 2 

Since the liquid-vapor interface is always at the saturation temperature, the average in- 
ternal energy of the liquid initially corresponds to the saturated temperature and pressure 
and at any later time is 


u 


i, av 





(D9) 


where 


m £ =m f,i- 4m v 


At any time the energy stored in the liquid must be equal to the sum of the energy stored 
in the two regions such that, 


m f u f , av " m b u b + m l 


% + u sat 


(DIO) 


The total mass of liquid is equal to the sum of the mass of liquid in the two regions. 


= m b + m l 


(Dll) 


Combination of equations (DIO) and (Dll) yields 


m b =m f 


^ u £, sat + “b ~ , ay ) 

u f, sat " u b 


(D12) 


All the liquid is initially saturated so the energy input to the bulk is 

Qb= m b,f u b,f- m i,i u l, S at,i (D13) 

Subtracting the energy input to the bulk from the total energy input to the liquid yields 


Q l " Q b 


(D14) 
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To determine the validity of the initial assumption, the percent filling of the inner sphere 
was recomputed 


pf= — X 100 (D15) 

^£,av^ 

where av is determined by using the computer and the values of the total sphere pres- 
sure and u £ ay as input. This new value of percent filling allows the vapor volume to be 
recomputed and the calculations are performed again starting with equation (Dl). This 
cycle is repeated until the liquid mass changes by less than 0. 1 percent; the computer 
then proceeds to the next time interval. 
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Figure 1. -Tank pressure as function of heat added for two 
energy distribution models. Tank size, 1 cubic foot 
(2.832x10-2 m3). 



Figure 2. -Liquid-hydrogen experiment apparatus. 






Figure 4. - Carbon resistor rakes. 


Temperature, K Pressure, P, N/cm : 


CVl 



(a) Total pressure as function of time. 



(b) Outer-sphere and average heater temperature 
as functions of time. 


Transducer 



(c) Upper inner-sphere internal and surface tem- 
perature as functions of time. 


(d) Lower inner-sphere internal temperature as 
function of time. 


Figure 5. - Test 1. Initial filling, 31.6 percent; average heat flux, 16.9 Btu per hour per square foot (53 W/rA (See fig. 3 
for location of transducers.) 
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Temperature, T, K Pressure, P, N/cnV 
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(a) Total pressure as function of time. 


(b) Outer-sphere and average heater temperature as 
functions of time. 
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Figure 6. - Test 2. Initial filling, 48.9 percent; average heat flux, 18.9 Btu per hour per square foot (60 W/m 2 ). (See 
fig. 3 for location of transducers. ) 
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(a) Total pressure as function of time. 
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(b) Outer-sphere and average heater temperature 
as function of time. 
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(c) Upper inner-sphere internal and surface tern- (d) Lower inner-sphere internal temperature as 

perature as functions of time. function of time. 

Figure 7. - Test 3. Initial filling, 79.8 percent; average heat flux, 21.9 Btu per hour per square foot (69 W/m 2 ). (See fig. 3 
for location of transducers. ) 
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(a) Total pressure as function of time. (b) Outer-sphere and average heater temperature 

as functions of time. 
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Figure 8. - Test 4. Initial filling, 54. 2 percent; average heat flux, 64.4 Btu per hour per square foot (203 W/m 2 ). (See 
fig. 3 for location of transducers.) 
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(a) 9-Inch (23-cm) diameter spherical tank (ref. 7). (b) 22-Inch (56-cm) diameter spherical tank. 

Figure 9. - Pressure as function of total heat added, effect of percent filling. 
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(a) 9-Inch (23 -cm) diameter spherical tank (ref. 7). 

Figure 10. - Pressure as 
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(b) 22-Inch (56-cm) diameter spherical tank. 
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Figure 11. - Inner sphere temperature as function of position for 

typical 22-inch (56-cm) diameter spherical tank test. 

Inner- 

sphere filling, approximately 50 percent. 
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Figure 12. - Pressure as function of total heat added divided by tank 
volume. (9-in. (23-cm) diam sphere data from ref. 7). 
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Figure 13. - Inner-sphere and heater emissivity as function of 
temperature. 



Figure 14. - Stainless-steel specific heat as function 
of temperature. (Data from ref. 20). 
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Figure 15, - Stainless-steel thermal conductivity 
as function of temperature. 
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